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It is well known that a small error in the measurement of the elastic constants will result, for all 
physically interesting boundary value problems, in small errors in the computed values of the stresses 
and displacements. In this paper actual bounds are given for the error in both the first and second 
boundary value problems. In addition it is shown that as Poisson's ratio tends to V2 the results for 
compressible theory tend to those for the analogous problems in the classical incompressible theory. 

1 . Introduction 

In the boundary value problems of classical elasticity the elastic constants of the material in 
question are assumed known. The solution of the problem is then expressed in terms of the 
"known" elastic constants. However, the values of these elastic constants are actually deter- 
mined empirically and are therefore not known exactly. Rather the values of these constants 
may be determined to lie between certain limits. 

It appears to be well understood by workers in classical elasticity that a small error in measure- 
ment of the elastic constants will not have much effect on the solution to the problem in question, 
and in problems where the solution can be determined explicitly for arbitrary elastic constants the 
effect of a small error can be easily determined. However, if the solution cannot be computed but 
instead pointwise bounds for the solution and its derivatives are sought, one must also take into 
account the error which might arise 1 due to inaccuracies in the determination of the elastic constants. 

In this paper we develop a priori inequalities appropriate to the first and second boundary value 
problems. These inequalities give upper and lower bounds for the pointwise error in approxima- 
tion of the exact solution by an arbitrary function. The error in measurement of the elastic 
constants as well as the error in approximation of the data of the problem is taken into account in 
the bounds. Thus by choosing trial functions to closely approximate the data, a close approxima- 
tion can be obtained, with a known bound on the error, provided the error in the constants is also 
known. The part of the error which results from the fact that the trial function does not satisfy 
exactly the data of the problem is quadratic. Thus, the Rayleigh-Ritz technique may be employed. 
This means that we choose from a class of trial functions the best function in the sense that the error 
term involving the data is minimized (see e.g., [3] 1 ). 

The above mentioned inequalities also are interesting in that they clearly display the fact that 
the solution to the first (or second) boundary value problem for an incompressible medium (Pois- 
son's ratio, cr, having the value V2) as formulated in the literature, is actually the limit of the solution 
for values of cr < V2. We show in fact that at points in a given region the solution for cr < V2 con- 
verges with all its derivatives to the classical solution. To the authors' knowledge this fact has 
not been pointed out in the literature. 



*lnstitute for Fluid Dynamics and Applied Mathematics, University of Maryland, College Park, Md.; consultants at the National Bureau of Standards, Wash- 
ington, D.C. 

1 Figures in brackets indicate the literature references at the end of this paper. 
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2. Notation 

Let D be a bounded domain with boundary T in three dimensions. The displacement vector 
u a (u^ u%, u%) corresponding to elastic constants /x a and <r a is assumed to satisfy in D the system 
of equations 

Au^pv-iuL-'Fi, £ = 1,2,3 

u j,j=- c «P a - (2- 1 ) 

Here A denotes the Laplace operator, the symbol, z", denotes partial differentiation with respect to 
Xi and the summation convention is employed for Latin but not Greek indices. The constant c a 
is given by 

c«=l-2o-«, (2.2) 

and Fj denotes the body force per unit volume. In terms of uf and the elastic constants the stress 
components rg are expressed as 

Tf j = ) Lta[<j+</-(l-Ca)p a 6o]. (2.3) 

The physically interesting values of c a lie in the interval 

0^c a <3. (2.4) 

In this paper we shall, for simplicity, restrict our attention to the first and the second boundary 
value problems in classical elasticity. In particular, we assume that in D 

Au«- P «. = -^Fr, ^-p^=-^Fi 

u lj = - c «P a ' ufj = -c^. (2.5) 

We wish to compare the solutions uf and uf if one of the following sets of boundary conditions is 
satisfied: 

u<* = u$=fi on L (2.5A) 

^~o X T?} n j = ^ l Tfj n j = gi on r (2.5B) 

Tfjiij = TfjJij = h; on T. (2.5C) 



From the point of view of actual physical applicability, problem (2.5A) and (2.5C) are the two 
interesting cases. 

In the subsequent sections we shall set 

Wi=u°l — uf\. (2.6) 



3. Point wise Convergence 

We demonstrate in this section that if p^ and Ft are square integrable in D and hi is square 
integrable on T, then at interior points of D, u™ and all its derivatives tend to u/f and its correspond- 
ing derivatives as c«— > c#, and pb a — > /x#. 
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We assume that there exists a constant Mp such that 

\jjPydv<M & . ,3.1) 



L 



Jn 

Let us assume now that we were able to obtain for either of problems (2.5A) to (2.5C), an 
L 2 bound for wi in terms of L 2 integrals of the data of uf. Suppose further that the coefficients of 
the integrals of data tend to zero as c a ~~ * cp and /ju a —+ /x/?. Then if the data of uf are in L> it follows 
that as a-^/3, Wi—* in L 2 . Pointwise convergence would then follow from the continuity of w\. 
The required L 2 bounds for w\ are derived in appendix 1 and appendix 2. Since we are also inter- 
ested in computing pointwise bounds, explicit constants will be obtained. In problem (2.5A) we 
have 

[ WiWidv ^ -^- (c« - c/0 2 + ^« 1 ~^fi 1)2 f FiFidv, (3 2) 

d 2X,c« A? % 

and in problem (2.5B) 

f tiwdt>*yF(c«)^ (3.3) 

These inequalities follow from (3.1) using (5.23) and (5.30) of appendix 1. For case (2.5C) we have 
instead from (3.1) and appendix 2, 

/,«"*« (2-3'm- (;-l) Ca] ^ + 4£ ^-ffH /r »'*** 

+ 8( g j L^s 1) *8f/«F«*. (3.4) 

where 7 satisfies (5.8). 

We wish to obtain from these expressions pointwise inequalities in terms of the appropriate 
sets of data. In order to do this we introduce the function Tf k (P, Q), the Kelvin fundamental solu- 
tion of (2.1); i.e., 

r?.(r,Q)= ] [ ( ] + 2 ^)^- 1 [xi(P)-x i (Q)][x,(P)-x k (Q)] 1 

lk ' 87r(l + Ca)[ r PQ rj* Q J 

where tpq is the distance from the point P to the point Q. Then a particular solution u? of (2.1) 
is given by 



u7(P) = fi? j D T?,(P, Q)F k dv Q . (3.6) 



Let uf now be defined as 



Clearly then ttf satisfies the system 



We now define 



uf = u?-uf. (3.7) 



A/1y + -#«,, = in D. (3.8) 



— A 

Wi = Wi + wt (3.9) 
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where 



7i = u<!-u(i = j D [MaTft(P,0) - tx?YUP,Q)Wkdv (3.10) 



and 

A 



Wi = u<?-&e. (3.11) 



A 



Since 3? and ftf are biharmonic functions, it follows that W\ is a biharmonic function. We may, 
therefore, make use of some mean value inequalities for biharmonic functions (see, e.g., Bramble 
and Payne [2]). That is, if R denotes the radius of the largest sphere contained in D with center 
at the point at which bounds are desired (we take this point to be the origin), then 

™<0)"«(0) ^Y^ 2 j D ™i™idv. (3.12) 

Similarly, as is shown in [3] , it is possible to obtain 

Z) v (w f )(0)D y A)(0) ^ C v J WiWidv. (3.13) 

where D v is any v\h order partial derivative and C v is a constant depending on R and v. 

It is clear at this point that if in any of problems (2.5A to 2.5C) the body force Ft is identically 

zero, then w = and w = wi. By using the appropriate bound for W[W\dv (given by (3.2), (3.3), 

Jd 

or (3.4) we would obtain the desired pointwise bounds for a? — uf. In fact, in any compact sub- 
domain of Z), we could obtain a uniform bound of the following type for any v\h order partial deriva- 
tive, D v (u^ — u^\ of uf — uf} (assuming that in (3.5), I hihtds is bounded): 

Z>-(u?-uf)^(^-af) l^o^^^Cce,-^] 2 ^^^^- 1 -^ 1 ] 2 . (3.14) 

In either of problems (2.5A) or (2.5B), the constant K v would in fact be zero. We turn now to the 
case Ft ^ 0. 
Clearly 



A A 



to l {0)w i (0) ^ 2[w i (0)w i (0) + ^(0)^(0)]. (3.15) 

But from (3.12) 



a a 75 



WiWidv~\- I 

JD J D 



WiWtdv+ WiWidv 



(3.16) 



8tt/? 3 
Combining (3.15) and (3.16) we obtain 

m(Q)wi(0) ^ t^ J wmdv + 2 ^(0)^(0) + -^ J wmdv. ( 3<17 ) 

We seek now an appropriate bound for the last two terms of (3.17). An application of Schwarz's 
inequality to (3.10) gives 



WiWi 



^t/^rs.-/^ ( 3.i8) 
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We now insert (3.5) into the first term on the right and make use of the arithmetic-geometric mean 
inequality to obtain 

( r «-r»-u-P'iru-r'-u-M>'w P <g^f2rM-'-M7'T + r (1+2ca) (1 + 2cg) TI 



where G(P)= rp\dvQ. 

Combining (3.18) and (3.19) and integrating we have 

r ^ i d^# i { 2 [^-^-T+[ ( \t! Ca \ - ,J 1 n! c A T] f ^**. 

Jd 327T 2 [ * L^a(l + c a ) /l^(1 + Cjb)J J Jd 

where // = I r^\dvqdvp. 
JdJd 

It follows now from (3.18), (3.19), and (3.20) that 



2wiiO)m(0) + ^|^ J w\ w,dv ■ 



1 Ag(0)+^t»h 



16tt 2 



2[/x- 1 "^ 1 ] 2 + 



for/? 3 

(l + 2c«) + 2 CjS ) 



JLt a (l-hc a ) ^(1 4-C/s) 

The following bounds on G(0) and // are easy to obtain: 



H 



G(0) ^ (4tt^) 



and 



(3.19) 



(3.20) 



F k F k dv. (3.21) 



(3.22) 



H^V(^jrd) 



(3.23) 



where F is the volume of D and J is the diameter. Thus combining (3.17), and (3.21), (3.22), (3.23), 
and (3.26) we obtain 



umm ^l w ^ d ^[t{ l+ ^ v )\ 



4t7# 



(l + 2cq) (1 + 2^) 



Ha(l + C a ) /JL a (l + C a ) 



F k F k dv. (3.24) 



If Fa- is square integrable then the insertion of (3.2), (3.3), or (3.5) into (3.24) gives in any compact 
subdomain D\ of D a uniform bound of the following type. 



[uf(P)-u^P)][u'!(P)-ui(P)]^B l [^- t i^y + B,[c a -c^ 



(3.25) 



where B\ and B> are independent of the point P in D, and can be chosen independent of the elastic 
constants. This follows because of the fact that since c a , cp 5= 0, 

+ 2c„) (1 



\ q+2c„) (l + 2cfl) 1 2 _r _,_ _ l+2c a _, (cg-Cff) 

L/i„(l + Ca) ft* 1 + C)\ i ,M ° M " 1 + C a I 1 * ( 1 + C„)( 1 + C, 



•e) 



2 {MtJi-'-^f + ^iCa-cnf}. (3.26) 



The inequality (3.25) clearly implies the uniform convergence of uf to uf in any compact subdomain 
of D as fx a —* up and c a ~ * eg. 
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It is clear that if Fk has a sufficient number of square integrable derivatives then the conver- 
gence of the derivatives of uf may be obtained in a similar way. To see this we note that 

drf k (p,Q) = drf k (P,Q) 

dxp 8xq 

Thus (3.6) may be differentiated with respect to the variables of the point P, the order of differentia- 
tion and integration being interchanged on the right. Then (3.27) is used followed by an integration 
by parts which takes the differentiation off of rf fc . The bounds for \uf [ P (P) | are then obtained in a 
manner similar to that used in obtaining (3.25). To obtain the desired expressions for higher 
derivatives, we successively differentiate Ff k and use (3.27) to throw the differentiation off of Ff k in 
the volume integral on the right. Then if a sufficient number of derivatives of Fa- are square 
integrable over D and over T we obtain a uniform bound similar to (3.25) for any desired derivative of 
u". This shows, then, that u" and its derivatives converge to u% and its derivatives. 

We note, in particular, that if cp = the uf problem is the classical boundary value problem 
for an incompressible medium. It follows then that as c a ^0 and /x a — >/x/3 the solution to the 
problem corresponding to a converges to the solution in the incompressible case. This fact, to the 
knowledge of the authors, has never been pointed out in the literature. 

Another interesting fact is that Fa need only have the desired number of square integrable 
derivatives in a neighborhood of the point at which convergence is sought. The procedures used in 
this section could be applied to any neighborhood ft of the origin. Then since 



I WiWidv ^ WiWjdv 
Jn Jd 



it would follow that if Fa were square integrable over D and enough, derivatives of Fa were square 
integrable over (1 and its boundary, uniform convergence in fi of the derivative would follow. 



4. Effect of Error in Measurement of Elastic Constants on Pointwise Bounds 

We consider now problem (2.1) with either the displacement components or the surface 
tractions prescribed on T. Since c a and fju a must be determined empirically we shall assume that 
some error was made in making the measurements and that the values c# and fip were actually 
obtained. If <Pi is a sufficiently smooth vector function chosen to approximate the solution u - we 
may write 

\u1-<pt\ ^|tt?-a?|H>f-<pi|. (4.1) 

We wish to obtain upper and lower bounds for uf(0). If we can make each term on the right of 
(4.1) less than some prescribed e we would then have 

-2e<uf-ipi<2e (4.2) 

and <pi would be a good pointwise approximation to uf. 

The bound for the first term on the right of (4.1) is obtained for the first boundary value prob- 
lem in elasticity from (3.24) and (5.22) of appendix 1. 

«(0)-«O))(T(O)-^(O))«^{^J / ^ + 3 " , -7 | ^y» }. 14.3, 

where 

^„/ i =(l/X? + 2/?'V/15 + 5^/47r)(At- | -^ 1 ) 2 + ^'(2/? :t /75+F/477-)(c a -c /J ) 2 . 
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A bound for E(uP, uP) is obtained by the methods of Diaz and Greenberg [4], Fichera [6], Synge 
[10], or Bramble and Payne [3, eq (3.12)]. 

The second term on the right of (4.1) may be dealt with by the method of Bramble and Payne 
[3, eq(4.10)]. 

We look now at the stress components rg. Using the triangle inequality we obtain 

kS~ Tij| « kg — rg| + |rg— tijI, (4.4) 

where 



Tij = ^n 






(4.5) 



is an approximation to rg. Bounds for the first term on the right of (4.4) may be obtained using 
(3.13) following again the method outlined at the end of section 3. The second term on the right 
may be dealt with as described in [3] . 

For the second boundary value problem in elasticity we again use (4.1) and (4.4), where in this 
case uf and uf are solutions to (2.5C) and bounds for the second term on the right are given in [1] . 

5. Appendix 1 

In appendix 1 we derive the inequalities (5.22), (5.23), (5.29), and (5.30) which are needed in the 
text. In order to do this we start with 

I Wi,j(wi,j + Wj,i)dv. 
Ji) 

In particular, we bound this expression in terms of the solution of the problem with elastic constants 
jLt/a, cp. Consider then the Green's identity 

Wi, j ( Wi, j + Wj, i )dv= I W\ ( Wij + Wj, i ) rijds — I w% ( Aw, + w h y x )ds. (5.1) 

Making use of (2.5) and integrating by parts we obtain for w x satisfying (2.6) 
I Wij{wi,j + Wj,i)dv=— I w i {n- ] T ( * j n j -^ l Tf j nj}ds-\- (u~ 1 -^ 1 ) j w { F\dv 

+ (c a +C{i-2cc£fi) I popddv-cail-ca) \ (p a ) 2 dv - Cfi(l - c^ \ (p*) 2 dv. (5.2) 
Jn Jd Ji) 

In either of problems (2.5A) or (2.5B), the boundary integral on the right of (5.2) will vanish. In 
fact, this integral would vanish for various mixed boundary conditions which are encountered in 
physical situations. On the other hand, in problem (2.5C) the boundary term does not vanish. We 
shall thus consider first problems (2.5A) and (2.5B) and then treat (2.5C) at the end of the section. 
If < c a ^ 3/4 (a ^ 1/8) we rewrite (5.2) for (A) and (B) as 



Wj,j(Wi,j + Wj, \)dv=— C a (l— Co) 
JD JI) 



(Cq + C/j — 2CqC/j) 



2 

dv 



2c a (l — Co) 

+ (Ma 1 -^ 1 ) f WiFidv+ £° { ~^ 2 a) J D (P*) 2 dv. (5.3) 
The first term on the right is nonpositive. Thus 

f wi,j(wij + Wj, i)dv ^ i C9 ( 7 C ^\ { (pvydv+ifi-'-^') f wiFidv. (5.4) 

JD *CaVl —C a ) JD JD 
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The strain energy E(u^, uP) corresponding to u 1 * is given by: 

2il?E(uP,uP) = [ [ufjiui^uld^c^l-c^ip^^dv^^ {^-c^) ^{jfiYdv. (5.5) 

The last inequality results from the fact that 

utjiufj + uf,,) «2[« I ) 2 +(uf, 2 ) 2 +(uf, 3 ) 2 ] ^2/3[<|. (5.6) 

Hence if c# ^ we have from (5.4) and (5.5) 

r 3lm~ x (c — cr) 2 f 

Wij(wij + Wj,i)dv^n 71 ^T7^ r^W, uP) + (p^ 1 — M^ 1 ) WjF } dv. (5.7) 

J/) ■ 2c a C^(l— C a ) (3 — Cfl) Jd 

Bounds for E(uP, u&) in problems (2.5A) and (2.5B) in terms of the data of the problem have been 
given for instance by Diaz and Greenberg [4], Fichera [6] , Prager and Synge [9], Synge [10], and 
Bramble and Payne [1,3]. If F\ = such bounds would give then for < c a ^ 3/4 an upper bound 
for 



\wi,j{wuj 



+ Wj, i)dv 



in terms of the data of the uf problem of (2.5A) or (2.5B). 
For c a > 3/4 we choose a constant y in the interval 

l-c" 1 <r<2/3 (5.8) 

and consider the following identity derived from (5.2) by addition of the appropriate term to each 
side of the equation 

Wi, ; ( Wi, } + Wj t i)dv — y wf dv= [c a + cp + 2 ( y — 1 ) c a ca\ p a p^dv 
Jd Jd Jd 

-c a {l + (y-l)c a } jjp«ydv-c«{l-(y-l)c«} j (p») 2 dv+ (p-* -p?) f wiF,dv. (5.9) 
Using (5.6) then we have 

(i-|)/^u(», J ^,^«-c. ! i + (,-i)^/Jp- + ^ 1 ^;_-'; ( ffV }^ 

+ 4^T^TM Jy ) '* + ^'-n')j^>' l A,. (5.101 

Thus we obtain 

The choice y = 2c a /(3 + 2c«) yields 

f w, M, i + wj, ddv ^ (3 \ 2c " )2(Ca ~ c ^ fiptydv + in- 1 -^ 1 ) fmFidv (5.12) 

JD ^Caio — Ca) JD P JD 

Proceeding as before we obtain finally the inequality 

Wij(wij + Wi,i)dv^F(c a )(Ca — cp) 2 (jfifdv + 0*" 1 ~ Mtf 1 ) Fmdv, (5.13) 

Jd Jd Jd 
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where 



I [4c«(l-c„)]- , ,0<c o «3/4 
F(c a ) = \ (5.14) 

[[(3 + 2c„) 2 [4c a (3-c a ) 2 ]- 1 ,c a >3/4. 

It follows then that for cp¥" 

f ^^(^j + ^j, i)rft; ^ ^'^^ ( Co -e^) 2 g(^, ^) -H (^Z 1 - /^^) f WiFidv. (5.15) 

Jd Cp(3 — Cfi) JD 

In problem (2.5A) we have the relation 

Wi,jWi,jdv = —\ WiAwidv = — I w i (p a — p (i )dv 

=- f (capo-cppP) (p a -pf*)dv+ (fJL- l -Pv l ) JFiWidv. (5.16) 
or 

f D Wi.jWi,jdv = - C a j D P*- ^*^ P^dv+itla 1 -*?) j^iWidv+ ^'^' j^ ( P ^dv.(S.l7) 

Thus, dropping the first term on the right, 

f w UJ w,,jdv^ ( Ca ~ Cp)2 f (pP) 2 dv+(tL- l -rf) ( Fmdv (5.18) 

Jd %c a Jd Jd 

and if c$ # 

f m.jm./ifa ^ 3^* (^~^) 2 g (^ ^) + (tf-rf) [ Fniwfc. (5.19) 

Jd zc a C/3(o — C&) J D 

We have retained the term (/x" 1 — /x^ 1 ) I FiWidv up to this point since in the important case 

Fi = this term will vanish and the inequality will not be altered. Let us assume now that Fi # 
and note that 

wmdv^—\ Wjjwijdv^ e — - — -& r— + — ^ ^— Wf- (5.20) 

Jd Xi Jo Zkic a cp(3 — c$) Xi Jd w.^w 

Where Xi is the first eigenvalue in the three dimensional fixed membrane problem for D. A lower 
bound for Xi is given by the Faber-Krahn inequality [5, 7], i.e., 

X,^[jj] 2/ \ (5.21) 

where V denotes the volume of D. From the arithmetic-geometric mean inequality applied to the 
last term of (5.20) it follows that 

f wmdv ^ \ (Ca ~% y *\ E( m i^)+ (M " 1 "^ 1)2 f FiFidv. (5 22) 

Jd XiCaCj3(3 — Cfi) Kf J D ^' l ^ 

Similarly if we insert (5.18) into the first inequality in (5.20) and apply the arithmetic-geometric 
mean inequality to the last term we obtain 

f wmdv « %=^1! f (p d)« A , + ^'-/ t ii')' f FiFidv , (5.23) 

Jd 2XiC a Jd Xf Jd v y 
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For problem (3.5B) it was shown in [1, eq (3.14)] that if w\ is normalized so that I w;,<is = 0then 

WiWids^-rj- \ Wi,j(wij-\-Wi,i)dv (5.24) 

Jr Zk 3 jD 

where k\ and /c 3 are constants given in [1, eqs (2.9) and (2.19)], which are computable in terms of the 
geometry of D. 

We now make use of the identity 

5 I WjWjdv = x k Wi{2riiWk + rikWi)ds + 2 I x k w% (w,, k + Wk t i)dv + 2 \ x k WkW%, \dv. (5.25) 
Jd Jr Jd Jd 

An application of the arithmetic geometric mean inequality and (3.8) then yields (see e.g., [1, eqs 
(2.26) to (2.28)]) 

f , 6r M f . , 28r 2 , f 

WjWjdv ^ —z- WiWjds-r — w x , j\Wi,j + Wj, i)dv (5.26) 

jd o Jr ^5 Jd 

where Tm denotes the maximum distance from the origin to T. From (5.24) and (5.26) it follows 
then that 

I WjWidv ^ k 4 I Wij(u)jj + Wj, i)dv (5.27) 

JD Jd 



where 

h-- 

Thus if in (2.5B) ca ^ 0, 



— rukilkn-}- ^F f\t 



(5.28) 



f Wiwidv ^ 12k ^^ Ca *£. E(uP, >) + (fi- 1 -^ 1 )^ f *7Ff<fo. (5.29) 

J/J C/j(^ — C/3j J/; 

Finally, after inserting (5.13) into (5.27) and using the arithmetic-geometric mean inequality on 
the last term we have 

I w i Widv^2k 4 F(ca){c a -c [ 3) 2 I (pPydv+k^bi- 1 — /x^ 1 ) 2 FjFjdv. (5.30) 

6. Appendix 2 

We derive here inequalities (6.3) and (6.4) which are needed to treat problem (2.5C). 

Again we start with Wij(wij + Wj, i)dv and instead of (5.2) we obtain 
jd 

I Wij(wi,j-\-Wj,i)dv = — \ Wihi[/JL- 1 — n^ds-^ (c a + C(] — 2caC(]) \ p a p$dv 

-C a (l-C a ) j^ ( P «) 2 <fo-C0(l-C0) J^ { P ^dv+ f^ 1 -^) £ «*F|<fo, (6.1) 

which follows from a calculation involving the use of the divergence theorem. We thus obtain as 
before 

I Wi tj (wuj + Wj,i)dv+ (2-3-,) [^a 1 -^ 1 ] j w t A,<is 

- (2^) ( "«-"''» J[ f ** ' 2(2- Wif + t'- Dc] 1 <*>» M 
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where y satisfies (5.8). It follows then from the arithmetic-geometric mean inequality, (5.24), and 
(5.27) that if Wi is normalized so that I Wids = then 



2 2. 



Wjj(wij + Wj,i)dv^-ri% — p v (ki/k 3 ) 2 hihids 
Jd 4>(Z — 5y) £ a Jr 



(Ca-Cfi) 2 f {p^V 
Jl> 



- + - 



2 (^-^ 



2(2-3y)[l+(y-l)c„] (2-3y) 
for any positive numbers a and b. In particular for a=%, 6 = x k (6.2) becomes 

I hjhjds 
Jr 



^-**//» 



(6.3) 
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As before inserting (6.3) into (5.27) we have 
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